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We study the exchange interactions in half-metallic Heusler alloys using first-principles calcula- 
tions in conjunction with the frozen-magnon approximation. The Curie temperature is estimated 
within both mean-field (MF) and random-phase-approximation (RPA) approaches. For the half- 
Heusler alloys NiMnSb and CoMnSb the dominant interaction is between the nearest Mn atoms. In 
this case the MF and RPA estimations differ strongly. The RPA approach provides better agree- 
ment with experiment. The exchange interactions are more complex in the case of full-Heusler alloys 
Co2MnSi and Co2CrAl where the dominant effects are the inter-sublattice interactions between the 
Mn(Cr) and Co atoms and between Co atoms at different sublattices. For these compounds we find 
that both MF and RPA give very close values of the Curie temperature slightly underestimating 
experimental quantities. We study the influence of the lattice compression on the magnetic proper- 
ties. The temperature dependence of the magnetization is calculated using the RPA method within 
both quantum mechanical and classical approaches. 

PACS numbers: 75.50.Cc, 75.30.Et, 75.30.Ds, 75.60.-d 



I. INTRODUCTION 

During the last decade the half- metallic ferromagnets 
have become one of the most studied classes of materials. 
The existence of a gap in the minority-spin band struc- 
ture leads to 100% spin-polarization of the electron states 
at the Fermi level and makes these systems attractive 
for applications in the emerging field of spintronicsi In 
half-metals the creation of a fully spin-polarized current 
should be possible that should maximize the efficiency of 
magnetoelectronics devices^ 

The half-metallicity was first predicted by de Groot 
and collaborators in 1983 when studying the band struc- 
ture of a half-Heusler alloy NiMnSb.- They found that 
the spin-down channel is semiconducting. In 2002 
Galanakis et al. have shown that the gap arises from the 
interaction between the d-orbitals of Ni and Mn creat- 
ing bonding and antibonding states separated by a gap. 4 
Ishida and collaborators have proposed that also the full- 
Heusler compounds of the type Co2MnZ, where Z stands 
for Si and Ge, are half-metals^ In these compounds the 
origin of half-metallicity is more complex than in the 
half- Heusler alloys because of the presence of the states 
located entirely at the Co sites. 6 Several other Heusler 
alloys have been predicted to be half-metals: 7 Akinaga 
and collaborators^ were able to crystallize a CrAs thin 
film in the zinc-blende structure, that is similar to the 
lattice of the Heusler alloys. The magnetic moment per 
formula unit was found to be close to 3/ib that corre- 
sponds to the integer value characteristic for half-metals. 
A number of further half-metallic materials are Cr02 in 
a metastable cubic phase, Fe304, the manganites (e.g 
Lao.7Sr .3Mn0 3 )S, the diluted magnetic semiconductors 
(e.g. Mn impurities in Si or GaAs)ii2*ii 



Besides strong spin polarization of the charge carriers 
in the ground state the spintronics materials must pos- 
sess a high Curie temperature to allow the applications 
in the devices operating at room temperature. Available 
experimental information shows that the Heusler alloys 
are promising systems also in this respect^ Up to now 
the main body of the theoretical studies was devoted 
to the properties of the half-metallic gap^ Recently, 
Chioncel and collaborators studied the influence of the 
correlation effects on the electron structure of CrAs^ 
They found that the spin-magnon interaction leads to 
the appearance of non-quasiparticle states in the spin- 
minority channel. The states are shown to lie above 
the Fermi level and to be sensitive to the value of the 
lattice constant. For a number of Heusler alloys it was 
shown that half-metallicity is preserved under tetragonal- 
ization of the crystal lattice^ and application of the hy- 
drostatic pressure^. Mavropoulos et al. studied the in- 
fluence of the spin-orbit coupling on the spin-polarization 
at the Fermi level and found the effect to be very small 17 
that is in agreement with a small orbital moment cal- 
culated by Galanakis^ Larson et <d~- have shown that 
the structure of Heusler alloys is stable with respect to 
the interchange of atoms and Orgassa and collaborators 
and Picozzi and collaborators have demonstrated that 
a small degree of disorder does not destroy the half- 
metallic gapi22i2i Dowben and Skomski have shown that 
at non-zero temperatures the spin-wave excitations lead 
to the presence at the Fermi level of the electron states 
with opposite spin projections leading to decreasing spin- 
polarization of the charge carriers^ 

Despite very strong interest to the half-metallic ferro- 
magnetism in Heusler alloys the number of theoretical 
studies of exchange interactions and Curie temperature 
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in Heusler alloys is still very small. The first contribution 
to the density functional theory of the exchange interac- 
tions in these systems was made in an early paper by 
Kubler et al>^ where the microscopic mechanisms of the 
magnetism of Heusler alloys were discussed on the basis 
of the comparison of the ferromagnetic and antiferromag- 
netic configurations of the Mn moments. Recently, the 
studies of the inter-atomic exchange interactions in sev- 
eral Heusler compounds were reported by the present au- 
thors and Kurtulus et q/i 24 i 25 i 26 . §a§ioglu et al. studied 
the exchange interactions in non- half- metallic Ni 2 MnZ 
(Z=Ga,In,Sn,Sb) and half-metallic Mn 2 VZ (Z=Al,Ge). 
The importance of the inter-sublattice exchange interac- 
tion has been demonstrated. For example, in the case 
of Mn 2 VZ (Z=Al,Ge) it was shown that the antiferro- 
magnetic coupling between the V and Mn moments sta- 
bilizes the ferromagnetic alignment of the Mn moments. 
Kubler^ estimated T c of NiMnSb to be 601 K to 701 
K depending on the approach used in the calculations. 
These values are in good correlation with experimental 
value of 730 K.±2 

The main task of the present contribution is the study 
of the exchange interactions in both half- and full- Heusler 
alloys. We use the calculated exchange parameters to es- 
timate the Curie temperature in both the random phase 
(Tq PA ) and the mean field approximations (T^ rFA ). In 
Section ^ we briefly discuss the formalism employed in 
the calculations. In Section II I II we present the results 
on the spin magnetic moments and the density of states 
(DOS) for four compounds studied: NiMnSb, CoMnSb, 
Co2MnSi and Co2CrAl. In Sect ion HVI we discuss the cal- 
culated exchange interactions and Curie temperatures. 
Section [V] is devoted to the consideration of the temper- 
ature dependence of magnetization. The films of Heusler 
alloys grown on different substrates can have different 
lattice parameters and, as a result, noticeable variation 
of the electron structure. Section [V] contains the sum- 
mary. In the Appendix, we present the formalism for the 
calculation of the Curie temperature of a multi-sublattice 
ferromagnet within the framework of the random phase 
approximation. 



II. CALCULATIONAL METHOD 

Half- and full-Heusler alloys crystallize in the Clb and 
L2i structures respectively (see Fig. The lattice con- 
sists from 4 interpenetrating fee lattices. In the case of 
the half-Heusler alloys (XYZ) one of the four sublattices 
is vacant. The Bravais lattice is in both cases fee. In full- 
Heusler alloy the atomic basis consists of four atoms. For 
example, in Co2MnSi the positions of the basis atoms in 
Wyckoff coordinates are the following: Co atoms at (000) 
and (i \ §), Mn at {\ \ ±), Si at (f § §). The Co atoms 
at the two different sublattices have the same local envi- 
ronment rotated by 90° with respect to the [001] axis. In 
half-Heusler compounds the position (| | |) is vacant. 

The calculations are carried out with the aug- 
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FIG. 1: Clf, and L2i structures adapted by the half- and full- 
Heusler alloys. The lattice consists from 4 interpenetrating fee 
lattices. In the case of the half-Heusler alloys (XYZ) one of 
the four sublattices is vacant. If all atoms were identical, the 
crystal structure would be a simple bec lattice 



mented spherical waves method (ASW) 28 within the 
atomic-sphere approximation (ASA)S The exchange- 
correlation potential is chosen in the generalized gradient 
approximation^ A dense Brillouin zone (BZ) sampling 
30 x 30 x 30 is used. The radii of all atomic spheres are 
chosen equal. In the case of half-Heusler alloys we intro- 
duce an empty sphere located at the unoccupied site. 



A. Exchange parameters 

The method for the calculation of exchange constants 
has been presented elsewhereiSi Here we give a brief 
overview. 

We describe the interatomic exchange interactions in 
terms of the classical Heisenberg Hamiltonian 



H eff = 



E 



E 

R,R' 
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In Eq. Q), the indices /j, and v number different sublat- 
tices and R and R' are the lattice vectors specifying the 
atoms within sublattices, is the unit vector pointing 
in the direction of the magnetic moment at site (/1, R). 

We employ the frozen-magnon approach to calculate 
interatomic Heisenberg exchange parameters^ The cal- 
culations involve few steps. In the first step, the exchange 
parameters between the atoms of a given sublattice \x are 
computed. The calculation is based on the evaluation of 
the energy of the frozen-magnon configurations defined 
by the following atomic polar and azimuthal angles 



y R 



^R 



q R 
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The constant phase § v is chosen equal to zero. The mag- 
netic moments of all other sublattices are kept parallel to 
the z axis. Within the Heisenberg model the energy 
of such configuration takes the form 



',qj 



sin 2 0JW(q) 
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where Eq m does not depend on q and the Fourier trans- 
form ./^(q) is defined by 

J^(q)=^J ^exp( iq .R). (4) 

R 

In the case of v = /i the sum in Eq. |@J does not 
include R = 0. Calculating E^^{9,c\) for a regular q- 
mesh in the Brillouin zone of the crystal and performing 
back Fourier transformation one gets exchange parame- 
ters Jq£ for sublattice \i. The determination of the ex- 
change interactions between the atoms of two different 
sublattices (J, and v is discussed in Ref. l24l 

B. Curie temperature 

The Curie temperature is estimated within two dif- 
ferent approaches: the mean-field approximation (MFA) 
and random phase approximation (RPA). The MFA for 
a multi-sublatticc material requires solving the system of 
coupled equation a 24 i 32 

where {e v ) is the average z component of ej^ and Jq " = 
Sr ^or- Eq EJcan be represented in the form of eigen- 
value matrix-problem 

(0 - TI)E = (6) 

where Ml/ = -^-Jq" i I is a unit matrix and E is the 
vector of (e v ). The largest eigenvalue of matrix gives 
the value ofT^ FA M. 

A more consequent method for the study of the ther- 
modynamics of Heisenberg systems is provided by the 
RPA approachiiS^i The RPA technique is intensively 
used for studies of both single-sublattic o 35 ! 36 a nd multi- 
sublattico 37 ! 38 ! 39 ! 40 ! 41 ! 42 ! 43 ! 44 ! 45 systems. In the case that 
only the exchange interactions within one sublattice are 
important the Curie temperature within the RPA is given 
by the relation 3 ^ 

1 = 31y i (7) 

k B T* PA 2 N 2f J(0) - J(q) ' {l) 

We use the RPA approach to study the temperature de- 
pendence of the magnetization in the temperature inter- 
val from K to Tq- The RPA technique for a multi- 
sublattice system is briefly presented in Appendix. 

III. DOS AND MAGNETIC MOMENTS 

A. NiMnSb and CoMnSb 

In this section we report the calculation of DOS and 
magnetic moments at different lattice parameters for 



TABLE I: Calculated atom-resolved and total spin moments 
in hb for NiMnSb, CoMnSb, Co 2 CrAl and Co 2 MnSi. All 
compounds are half-metallic at the experimental lattice con- 
stants taken from Ref. llH an means the use of the lattice 
constant that places the Fermi level at the upper edge of the 
half-metallic gap and am corresponds to 1% contraction of 
the lattice constant with respect to a//. 



Compound 


a(A) 


X 


Y 


Z 


Void 


Total 


NiMnSb - 


&I[exp] 


5.93 


0.20 


3.85 


-0.09 


0.04 


4.00 


NiMnSb - 


an 


5.68 


0.32 


3.68 


-0.05 


0.05 


4.00 


NiMnSb - 


am 


5.62 


0.33 


3.64 


-0.04 


0.05 


3.97 


CoMnSb - 


&I[exp] 


5.87 


-0.32 


3.41 


-0.11 


0.02 


3.00 


CoMnSb - 




5.22 


0.45 


2.57 


-0.06 


0.04 


3.00 


CoMnSb - 


am 


5.17 


0.48 


2.52 


-0.05 


0.04 


2.99 


Co 2 CrAl - 




5.74 


0.62 


1.83 


-0.08 




3.00 


Co 2 CrAl - 


an 


5.55 


0.69 


1.68 


-0.06 




3.00 


Co 2 CrAl - 


am 


5.49 


0.69 


1.66 


-0.05 




2.99 


Co 2 MnSi - 


^I[exp] 


5.65 


0.93 


3.21 


-0.06 




5.00 


Co 2 MnSi - 


- an 


5.49 


0.97 


3.10 


-0.04 




5.00 


Co 2 MnSi - 


- "am 


5.43 


0.97 


3.01 


-0.04 




4.97 



NiMnSb and for CoMnSb compound that has one elec- 
tron per formula unit less than NiMnSb. The electronic 
structure of both compounds has been extensively stud- 
ied earlier and the reader is referred to the review 
and references therein for detailed discussion. Here we 
present a brief description of the calculational results 
aiming to provide the basis for further considerations and 
to allow the comparison with previous work. 

In Table[I]we collect the atomic and total spin moments 
for three different lattice parameters. The investigation 
of the influence of the value of the lattice parameter on 
the properties of the Heusler alloys is important since the 
samples grown on different substrates can have different 
lattice spacings. The first calculation is performed for 
the experimental bulk lattice constant M~ The calculated 
densities of states (DOS) for this case are presented in the 
upper panel of Fig. For both NiMnSb and CoMnSb 
the Fermi level lies in the low-energy part of the half- 
metallic gap. The compression of the lattice pushes the 
majority p states to higher energies that results in in- 
creased energy position of the Fermi level with respect 
to the half-metallic gap. At the lattice parameter an 
the Fermi level coincides with the upper edge of the gap 
(Fig. [21). In the next step we further contracted the lattice 
constant by 1% (lattice parameter am, bottom panel in 
Fig. |2J| . In this case the Fermi level is slightly above the 
gap and the total spin moment is slightly smaller than 
the integer values of 3 and 4 /ie for CoMnSb and NiMnSb 
respectively. 

The contraction of the lattice leads to an increase 
of the hybridization between the d orbitals of different 
transition-metal atoms. This results in a decrease of the 
spin moment of Mn. In the case of NiMnSb this change 
is small: the reduction of the Mn spin moment under 
lattice contraction from the experimental lattice param- 
eter to an is ^0.2 /is- The Ni spin moment increases by 
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FIG. 2: (a) Calculated spin-resolved density of states of NiMnSb and CoMnSb for three values of the lattice parameter, (b) 
The same for Co2CrAl and Co2MnSi. The upper panels present the results for the experimental lattice constant!-" The middle 
panels show the results for the lattice parameter a,u that is determined by the coincidence of the Fermi level with the upper edge 
of the half-metallic gap. The bottom panels present the results for lattice parameter a.m that is obtained by a 1% contraction 
of an . 



about the same value to preserve the integer value of the 
total spin moment of 4 /zb- 

In CoMnSb, the half-metallic gap is larger than in 
NiMnSb. As a result, the transition of the Fermi level 
to the upper gap-edge requires a large lattice contraction 
of 11% (Table P). This leads to a strong decrease of the 
Mn moment by 0.84^b- To compensate this decrease the 
Co moment changes its sign transforming the magnetic 
structure from ferrimagnetic to ferromagnetic. 

The influence of the lattice contraction on the exchange 
interactions and Curie temperature is discussed in the 
next Section. 



reflects a smaller value of the Cr moment compared to 
the Mn moment (Table D). 

As in the case of the half-Heusler compounds discussed 
above, the variation of the lattice parameter leads to the 
change in the position of the Fermi level. At the ex- 
perimental lattice parameter the Fermi level of Co2CrAl 
lies in the lower part of the half-metallic gap while for 
Co2MnSi it is close to the middle of the gap (Fig. I^Jd). 
The contraction of the lattice needed to place the Fermi 
level at the upper edge of the gap is smaller than for 
CoMnSb. As a result, the change in the magnetic mo- 
ments is also relatively weak (Table QJ. 



B. Co 2 CrAl and Co 2 MnSi 

The second group of materials studied in the paper 
is formed by the full-Heusler compounds Co2MnSi and 
Co2CrAl. The electronic structure of these systems has 
been studied earlier Compared to half-Heusler systems, 
the presence of two Co atoms per formula unit results 
in an increased coordination number of Co atoms sur- 
rounding Mn atoms (eight instead of four in CoMnSb). 
This leads to an increased hybridization between the 3d 
orbitals of the Mn and Co atoms. The spin moment of 
Co in Co2MnSi is about 1 /ib that is considerably larger 
than the Co moment in CoMnSb. In Co2CrAl the Co 
moment is about l/3rd smaller than in Co2MnSi that 



IV. EXCHANGE PARAMETERS AND CURIE 
TEMPERATURE 

A. NiMnSb and CoMnSb 

In Fig. [3]we present the exchange constants calculated 
for various lattice spacings. The Co-Co, Ni-Ni exchange 
interactions as well as the exchange interactions between 
the moments of the 3d atoms and the induced moments of 
Sb atoms are very weak and are not shown. The weakness 
of the effective Co-Co and Ni-Ni exchange interactions 
can be explained by a relatively large distance between 
atoms (Fig. ^) and relatively small atomic moments. 

On the other hand, each Ni(Co) atom is surrounded 



5 



by four Mn atoms as nearest neighbors that results in 
strong Mn-Ni(Co) exchange interaction (Fig. |3J). Also 
the exchange interaction between large Mn moments is 
strong. 

The ferromagnetic Mn-Mn interactions are mainly re- 
sponsible for the stable ferromagnetism of these materi- 
als. For both systems and for all lattice spacings stud- 
ied the leading Mn-Mn exchange interaction is strongly 
positive. In NiMnSb, the Mn-Ni interaction of the near- 
est neighbors is positive for all three lattice parameters 
leading to the parallel orientation of the spins of the Mn 
and Ni atoms. In CoMnSb the situation is different. At 
the experimental lattice parameter the leading Mn-Co 
interaction is negative resulting in the ferrimagnetism of 
the system. For the contracted lattices the interaction 
changes sign resulting in the ferromagnetic ground state 
of the alloy. 

The analysis of the strength of the exchange interac- 
tion as a function of the lattice parameter shows that 
in CoMnSb the contraction leads to a strong increase of 
both leading Mn-Co and Mn-Mn interactions. On the 
other hand, in NiMnSb the increase of the Mn-Ni in- 
teraction is accompanied by a decrease of the leading 
Mn-Mn interaction. Simultaneously, the interaction be- 
tween the second-nearest Mn atoms increases with con- 
traction in the case of NiMnSb staying almost unchanged 
in CoMnSb. This complexity of the behavior reflects the 
complexity of the electronic structure of the systems. 

The interatomic exchange parameters are used to eval- 
uate the Curie temperature within two different ap- 
proaches: MFA and RPA. In Table |n] we present the 
values of the Curie temperature obtained, first, by taking 
into account the Mn-Mn interactions only and, second, 
with account for both Mn-Mn and Mn-Ni(Co) interac- 
tions. The contribution of the intcr-sublattice interac- 
tions to the Curie temperature appears to be less than 
5 % for both compounds and the Curie temperature is 
mainly determined by the intra-sublattice Mn-Mn inter- 
action. 

The MFA and RPA estimations of the Curie tempera- 
ture differ rather strongly fTable lTT|) . The relative differ- 
ence of two estimations is about 20%. The reason behind 
this difference will be discussed in the following section. 
For the systems considered here the RPA estimations of 
the Curie temperatures are in good agreement with the 
experiment, somewhat overestimating the experimental 
values. 

Recently Kiibler 27 reported estimations of the Curie 
temperature of NiMnSb. His approach is based on the 
evaluation of the non-uniform magnetic susceptibility on 
the basis of the Landau-type expansion for the free en- 
ergy. Within some approximations the parameters used 
in the study of the thermodynamical properties can be 
expressed in terms of the quantities evaluated within the 
first-principles DFT calculations. The estimated values 
of the Curie temperature are 601 K for a static approach 
and 701 K if the frequency dependence of the suscepti- 
bility is taken into account. These estimations are some- 
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FIG. 3: The variation of the interatomic exchange parame- 
ters for NiMnSb (upper panel) and CoMnSb (bottom panel) 
as a function of the interatomic distance. The left panel cor- 
responds to the experimental lattice constant, the middle and 
right panels correspond respectively to an and am parame- 
ters. 



TABLE II: Calculated Curie temperatures. The second and 



third columns contain the T, 



MFA(RPA) 



obtained with the ac- 



count for Mn-Mn (Cr-Cr) interactions only. In the next two 
columns all interactions are taken into account. The last col- 
umn presents the experimental values of the Curie tempera- 
ture from Ref. IT^ 



Tc (K) 




MFA-Y 


RPA-Y 


MFA-all 


RPA-all 


Exp. 


NiMnSb - 


^I[exp] 


1096 


880 


1112 


900 


730 


NiMnSb - 


8-II 


1060 


853 


1107 


908 




NiMnSb - 


&iii 


1008 


802 


1063 


869 




CoMnSb - 


8j[exp] 


785 


619 


815 


671 


490 


CoMnSb - 




1185 


940 


1276 


1052 




CoMnSb - 


am 


1140 


893 


1252 


1032 




Co 2 CrAl - 


&I[exp] 


148 


141 


280 


270 


334 


Co 2 CrAl - 


&I I 


168 


159 


384 


365 




Co 2 CrAl - 


a,m 


164 


154 


400 


379 




Co 2 MnSi 


&I[exp] 


232 


196 


857 


740 


985 


Co 2 MnSi - 


8.1 1 


142 


118 


934 


804 




Co 2 MnSi - 


8m 


110 


75 


957 


817 





what lower than the value of 880 K given by the RPA 
approach ( Table ITTf) . A detailed comparative analysis of 
the two calculational schemes is needed to get an insight 
in the physical origin of this difference. 

The contraction of the lattice in the case of the NiMnSb 
compound leads to an increase of the Mn-Ni interactions 
(Fig. |3J) . This results in increased difference between the 
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FIG. 4: (a) The exchange constants for Co2CrAl as a function of the interatomic distance, (b) The same for Co2MnSi. The left 
panels correspond to the experimental lattice constant, the middle and right panels to an and am parameters respectively. 
The superscripts 1 and 2 denote Co atoms belonging to different sublattices (Fig. . For comparison, the exchange parameters 
of Co2MnSi obtained in Ref. [2(] at the experimental lattice parameter (an exp ]) are shown. The following symbols are used in 
the presentation: + for the Mn-Mn and Mn-Co interactions and * for the Co-Co and Co 1 -Co 2 interactions. 



Curie temperatures calculated with the Mn-Mn interac- 
tions only and with both Mn-Mn and Ni-Mn interactions 
taken into account (Table ITT1) . For CoMnSb, the leading 
exchange interactions of both Mn-Mn and Mn-Co types 
increase in the value under transition from the experi- 
mental lattice constant to ajj (Fig. |3J). As a result, the 
Curie temperature increases with contraction by about 
50%. 



B. Co 2 CrAl and Co 2 MnSi 

The presence of an extra Co atom in the full-Hcuslcr 
alloys makes the interactions more complex than in the 
case of the half-Heusler alloys. In CoMnSb the impor- 
tant interactions arise between nearest Mn atoms (Mn- 
Mn interactions) and between nearest Mn and Co atoms 
(Mn-Co interaction). In the case of Co 2 MnSi (Fig. 0} 
the interactions between Co atoms at the same sublat- 
tice (Co-Co) and between Co atoms at different sublat- 
tices (Co 1 -Co 2 ) must be taken into account. The cobalt 
atoms at different sublattices have the same local envi- 
ronment rotated by 90° about the [001] axis. The leading 
interaction responsible for the stability of the ferromag- 
netism is the Mn-Co interaction between Mn atoms and 
eight nearest Co atoms (Fig. . This interaction changes 
weakly with the contraction of the lattice. Our exchange 



parameters agree well with the parameters of Kurtulus 
et al. (Fig. 2} who also found the Co-Mn exchange inter- 
action to be leading^ 

The interaction between nearest Co atoms at different 
sublattices (empty squares in Fig.^J) favors the ferromag- 
netism also and is stronger than the ferromagnetic inter- 
action between the nearest Mn atoms (filled spheres) . Al- 
though the spin moment of Mn atoms is larger than the 
moment of Co atoms (Table HJ the opposite relation be- 
tween exchange parameters can be the consequence of the 
smaller distance between the Co atoms: a/2 between the 
Co atoms and \/2a/2 between the Mn atoms. An inter- 
esting feature of the intra-sublattice Mn-Mn and Co-Co 
interactions is different signs of the exchange parameters 
for different distances between atoms. This leads to a 
RKKY-like oscillations of the parameters (Fig. HJ. 

In Co2CrAl the leading Cr-Co interactions (filled tri- 
angles) are much smaller than corresponding Mn-Co in- 
teractions in Co2MnSi. On the other hand, the lead- 
ing inter-sublattice ferromagnetic Co-Co interactions are 
comparable in both systems. The compression of the 
lattice leads to an increase of the magnitude of the 
inter-sublattice Co-Cr and Co 1 -Co 2 coupling. The intra- 
sublattice Cr-Cr and Co-Co interactions oscillate with 
varying inter-atomic distances. 

The difference in the properties of the exchange pa- 
rameters of the half- and full-Heusler alloys is reflected 
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FIG. 5: The frozen-magnon dispersions for NiMnSb and 
Co2CrAl along the symmetry lines in the Brillouin zone. 
The magnons correspond to the Mn sublattice in the case 
of NiMnSb and to the Cr sublattice in the case of Co2CrAl. 
For Co2CrAl, the frozen-magnon energies are multiplied by a 
factor of 8. 
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in the calculated Curie temperatures ( Table ITT|l . In con- 
trast to CoMnSb where the Mn-Mn exchange interactions 
are dominant, in Co2MnSi they play a secondary role. 
The rY^FAiRPA) ca i cu i a ^ ec j taking into account these in- 
teractions only is much smaller than the Curie temper- 
ature calculated with all inter-atomic exchange interac- 
tions taken into account fTableHTjl. The same conclusion 
is valid for Co2CrAl where the Cr-Cr interactions give 
about half of the Curie temperature obtained with all 
interactions included into consideration. 

A striking feature of the full Heusler compound 
Co2CrAl that differs it strongly from the half- Heusler sys- 
tems considered in the previous Section is a very small 
difference between the Tc values calculated within the 
MFA and RPA approaches. A similar behavior was ob- 
tained for the Curie temperatures of the zincblende MnSi 
and MnC4& In Co2MnSi, the relative difference of the 
MFA and RPA estimations assumes an intermediate po- 
sition between the half-Heusler systems and Co2CrAl. 

To understand the origin of the strong variation of the 
relative difference of the MFA and RPA estimations of 
the Curie temperature we compare in Fig. [S] the frozen 
magnon dispersions for two compounds. The magnons 
correspond to the Mn sublattice in the case of NiMnSb 
and to the Cr sublattice in the case of Co2CrAl. As 
seen from Table [H] the MFA and RPA estimations ob- 
tained with the use of these dispersions differ by 20% for 
NiMnSb and by 5% for Co 2 CrAl. 

The Curie temperature is given by the average value 
of the magnon energies. In MFA this is the arithmetic 
average while in RPA this is harmonic average. There- 
fore we need to understand why for Co2CrAl these two 
averages are much closer than for NiMnSb. The follow- 
ing properties of the averages are important for us. The 



FIG. 6: The spin wave dispersions for NiMnSb and Co2CrAl 
along the symmetry lines in the Brillouin zone. Filled and 
empty spheres denote the acoustic and optical branches, re- 
spectively. 



arithmetic average takes all the magnon values with equal 
weight whereas in the harmonic average the weight de- 
creases with increasing energy of the magnon. 35 '— '—•—- 48 
It is an arithmetic property that the MFA estimation is 
larger than the RPA one or equal to it if all numbers to 
be averaged are equal to each other. In terms of magnon 
energies, T^f FA is equal to Tg PA in the case that the 
magnon spectrum is dispersion-less. 

Considering the frozen-magnon dispersions from the 
viewpoint of these properties we indeed can expect that 
the arithmetic and harmonic averages will be closer for 
Co2CrAl. In Fig.[S]both curves are scaled to have almost 
the same maximal value. It is seen that the Co2CrAl 
dispersion has smaller relative contribution of the low- 
energy magnons because of the steeper increase of the 
curve at small wave vectors. It has also smaller contri- 
bution of the magnons with the largest energies because 
the maxima have the form of well-defined peaks oppo- 
site to NiMnSb where we get a plateau. Thus the main 
contribution in the case of Co2CrAl comes from interme- 
diate energies that makes the MFA and RPA estimations 
closer. 

In Fig. HO we present the calculated spin-wave spec- 
tra for NiMnSb and Co2CrAl. The spin-wave energies 
are obtained by the diagonalization of the matrix of ex- 
change parameters that contains all important intra- and 
inter-sublattice interactions. The number of branches in 
the spectrum is equal to the number of magnetic atoms 
in the unit cell: two in NiMnSb and three in Co2CrAl. 
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One of the branches is acoustic and has zero energy for 
zero wave vector. Also in the spin-wave spectra, we see 
strong difference between two systems. In NiMnSb, the 
acoustic branch is predominantly of the Ni type stem- 
ming from the weak interaction between Mn and Ni mag- 
netic moments (see Fig. OJ. On the other hand, the 
optical branch is of predominantly the Mn type. The 
strong hybridization between two sublattices is obtained 
only about q = 0. In Co2CrAl, the energy scale of the 
branches differs much smaller and the hybridization be- 
tween sublattices is stronger than for NiMnSb. 

Coming back to the considerations of the Curie tem- 
peratures, we conclude that, in general, the Curie tem- 
peratures of Co2MnSi and Co2CrAl calculated within 
both MFA and RPA are in good agreement with exper- 
iment while the MFA values in the case of NiMnSb and 
CoMnSb overestimate the Curie temperature strongly. 

The lattice contraction leads in both compounds to an 
enhancement of the Mn-Co(Cr-Co) exchange constants 
that results in an increase of the Curie temperature. 

Kurtulus and collaborators have calculated the Curie 
temperature for Co2MnSi within MFA and found the 
value of 1251 K that is considerably larger than our MFA 
estimate of 857 K. This difference is unexpected since the 
values of the exchange parameters obtained by Kurtulus 
et al. agree well with our parameters (Fig^J. To re- 
veal the origin of the discrepancy we performed the MFA 
calculation of the Curie temperature with the exchange 
parameters of Kurtulus et al. and obtained the Tq value 
of 942 K which is in reasonable agreement with our es- 
timate. Apparently the reason for the inconsistency is 
in the procedure of the solving of the multiple-sublattice 
MFA problem used by Kurtulus et al. that should devi- 
ate from the standard one^ 



V. TEMPERATURE DEPENDENCE OF THE 
MAGNETIZATION 

The study of the temperature dependence of the mag- 
netic properties of itinerant ferromagnets is one of the 
fundamental problems of ongoing researches. Although 
density functional theory can formally be extended to 
the finite temperatures^, it is rarely used because of 
the lack of suitable exchange-correlation potentials for 
magnetic systems at finite temperatures. Statistical me- 
chanics treatment of model Hamiltonians is usually em- 
ployed. In this section we will present the results of the 
calculation of the temperature dependence of magnetiza- 
tion that is based on the consideration of the Heisenberg 
hamiltonian with exchange parameters calculated within 
a parameter-free DFT approach (Sect. ITT|) . 

To calculate the temperature dependence of the mag- 
netization we use the RPA method as described in ap- 
pendix. We consider both classical-spin and quantum- 
spin cases. 

In the classical-spin calculations the calculated values 
of the magnetic moments (TableHJ are used. To perform 



quantum-mechanical RPA calculation we assign integer 
values to the atomic moments. In the semi Heusler com- 
pounds we ignore the induced moments on Ni and Co 
atoms and assign the whole moment per formula unit 
to the Mn atom: 4/i B (S = 2) in NiMnSb and 3fi B 
(5 = 3/2) in CoMnSb. In Co 2 MnSi we take the values of 
3^3 (S = 3/2) and lfi B {S = 1/2) for Mn and Co atoms 
respectively. This assignment preserves the value of the 
total spin moment per chemical unit. In Co2CrAl we use 
in the quantum-RPA calculations the atomic moment of 
2fj, B (S = 1) for Cr and lfi B (S = 1/2) for Co. 

In Fig. Ufa), we present in the normalized form the cal- 
culated temperature dependence of the magnetization for 
both families of Heusler compounds. The calculations are 
performed for the experimental lattice parameter. For 
comparison, the experimental curves are presented. The 
nature of the spin (quantum or classical) influences the 
form of the curves considerably. The classical curve lies 
lower than the quantum one. This results from a faster 
drop of the magnetization in the low-temperature region 
in the case of classical spins. In general, the quantum 
consideration gives better agreement of the form of the 
temperature dependence of the magnetization with ex- 
periment. 

In Fig. [7[b) we present the temperature dependence of 
the magnetization of individual sublattices. As expected 
from the previous discussions in half-Heusler systems the 
main contribution to the magnetization comes from the 
Mn sublattice while for the full-Heusler systems both 3d 
atoms contribute substantially. 

Considering the calculated Curie temperatures we no- 
tice that the value of Tq calculated within the quantum- 
mechanical RPA is substantially larger than the corre- 
sponding classical estimation (see Fig. 0. This prop- 
erty is well-known and has its mathematical origin in 
the factor (S + 1)/S entering the RPA expression for the 
Curie temperature (Eq. IA14jl . In Fig. |S| we show the 
dependence of the Curie temperature calculated within 
the quantum mechanical RPA approach on the value 
of S. The exchange parameters are kept unchanged in 
these calculations. We see that the dependence has a 
monotonous character tending to a classical limit for 
large S. 

Presently we do not have an explanation why 
quantum-mechanical calculations give better form of the 
temperature dependence while the classical calculation 
provides better value of the Curie temperature. We can 
suggest the following arguments. The quantum treat- 
ment is more appropriate than the classical one in the 
low-temperature region. At high temperatures charac- 
terized by strong deviation of the atomic spins from the 
magnetization axis the quantum treatment gives too slow 
decrease of the magnetization. It is worth noting that the 
consequent theory should take into account not only the 
orientational disorder of the atomic moments but also 
the single-particle (Stoner-type) excitations leading to 
the decrease of atomic moments. Another important as- 
pect is related to the fact that the exchange parameters 
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FIG. 7: (a) The calculated temperature dependence of the total magnetization for both families of Heusler alloys. For com- 
parison the experimental temperature dependences— are presented. The calculations are performed for both classical and 
quantum Hamiltonians. Both the magnetization and the temperature are given in reduced form, (b) Calculated sublattice 
magnetizations as a function of temperature. The temperature is given in reduced form. 



used in the calculations are estimated within the picture 
of classical atomic moments described above. It is pos- 
sible that the values of the exchange parameters must 
be modified for the use in the quantum-mechanical case. 
These questions belong to fundamental problems of the 
quantum-mechanical description of the magnetic systems 
with itinerant electrons. 



VI. SUMMARY AND CONCLUSIONS 

We studied the electronic structure of several Heusler 
alloys using the augmented spherical waves method in 
conjunction with the generalized gradient approxima- 
tion to the exchange and correlation potential. Using 
the frozen-magnon approximation we calculated inter- 
atomic exchange parameters that were used to estimate 
the Curie temperature. The Curie temperature was esti- 
mated within both mean-field and random-phase approx- 
imation techniques. 

For the half-Heusler alloys NiMnSb and CoMnSb the 
dominant interaction is between the Mn atoms. The lat- 
tice compression results in considerable change of the ex- 
change parameters and Curie temperature. 

The magnetic interactions are more complex in full- 
Heusler alloys Co2MnSi and Co2CrAl. In both cases the 
ferromagnetism is stabilized by the inter-sublattice inter- 
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FIG. 8: Curie temperature of NiMnSb and CoMnSb as a 
function of spin quantum number S. The horizontal broken 
lines correspond to the classical limit (5* = oo). 



actions between the Mn(Cr) and Co atoms and between 
Co atoms belonging to different sublattices. Both the 
random phase and mean field approximations slightly un- 
derestimate the values of the Curie temperature. Com- 
pression of the lattice constant has little effect on the 
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magnetic properties of the full-Heusler alloys. 

We study the temperature dependence of the magneti- 
zation within the quantum mechanical and classical RPA. 
The quantum-mechanical approach gives the form of the 
temperature dependence that is in good agreement with 
experiment. The value of the Curie temperature is, how- 
ever, overestimated in the quantum-mechanical calcula- 
tion. 



Acknowledgments 

The financial support of Bundcsministerium fur Bil- 
dung und Forschung is acknowledged. I.G. is a fellow of 
the Greek State Scholarship Foundation. We thank the 
authors of Ref. for making the manuscript available 
before publication. 



APPENDIX A: THE RANDOM PHASE 
APPROXIMATION FOR MULTI-SUBLATTICE 
HEISENBERG HAMILTONIAN 

The Green function approach is a powerful tool 
in the study of the magnetism of complex sys- 
tems. (See, e.g., the application of the method to 
antifcrromagnet^l, ferrimagnets^&i^, random alloys 4 ^, 
layered system o 41 i 42 i 43 , disordered dilute magnetic 
systems^, multi-sublattice ferromagnets^.) In this ap- 
pendix we briefly overview the formalism to study the 
temperature dependence of the magnetization of multi- 
sublattice systems within the random phase approxima- 
tion. 

We start with the Heisenberg Hamiltonian for quan- 
tum spins 

ij p,v 

where e ijA( = (sf^, sf l( J/(S>) is the normalized spin 
operator corresponding to site 

In terms of the creation and destruction operators 
I? = sf =F s^„ the Hamiltonian can be written in the 
form 

ff =-EECte + « (A2) 

ij H" 

where j£? = Jf? /SpS v . 

Following Callen^i let us introduce Green function 

Gt;(r) = ~e(r){[sfJr),eMv§lu)sJ, v }) (A3) 

where 77 is a parameter, 9(t) is the step function (9(t) = 1 
for t > 0), [. . .] denotes the commutator and (. . .) is the 
thermal average over the canonical ensemble, ie., (F) — 
Tr[exp(-/3ff)F]/Tr[exp(-/?i0] with (3 = l/k B T 



Writing the equation of motion for G^(t) we obtain 

x{[[iW,ff],exp(t)sys:J) (A4) 

The last commutator term in Eq. I|A4[) generates higher- 
order Green functions. These functions can be reduced 
to lower-order functions by using Tyablikov decoupling 
(random phase approximation) scheme^: 

ra (^)([4Ar),sT v }) (A5) 

Applying this decoupling procedure to Eq. (IA4|) we get 

^ -^W<[<pW.exp(^,v)^J> 

+f EW^M 

-{sU)G^{r)\ (A6) 

After a Fourier transformation to energy and momentum 
space [g(q,uj) = j^J2iI dtue~ lclR 'Gio(T)] we obtain 

Tiugnv(<l,v) = 2~([*J,exp(77s^)s~])5 MV 

- 2 ^{^(q)(^).9^(q> w ) 

s 

-■^(0){^, c )s^(q,w)} (A7) 
Eq. I|A7() can be written in a compact matrix form 

M-M(q)]g(q,w)=u (A8) 

where g(q, ui) is a symmetric square matrix, I is a unit 
matrix and the inhomogeneity matrix u is expressed by 

*V = ^([^ ex P(^)C])V> (A9) 



matrix M(q) is defined by 




Next, we introduce the transformation which diagonal- 
izes matrix M(q)44 

L(q)M(q)R(q) = fi(q) (All) 

where f2(q) is the diagonal matrix whose elements give 
the spin wave energies a> M (q). The number of branches in 
the spin wave spectrum is equal to the number of mag- 
netic atoms in the unit cell. The transformation matrix 
R(q) and its inverse R _1 (q) = L(q) are obtained from 
the right eigenvectors of M(q) as columns and from the 
left eigenvectors as rows, respectively. 
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Using the spectral theorem and Callen's techniques 
one obtains the thermal averages of the sublattice mag- 
netizations: 

= (s, - ggXi + ^) 2S - +1 + (s, + 1 + <M<^' +1 

(1 + _ ($ p )2S M +l 

(A12) 

where $ M is an auxiliary function given by 

^ = ^EEV(q) e ^ ( 1 q) _ 1 ^(q) (Ai3) 

q v 

In Eq. I|A13(I . N is the number of q points in the first 
BZ. 

Eq. i|A12l) is the central equation for the calculation 
of the sublattice magnetizations. It must be solved self- 
consistently. The Curie temperature Tq is determined as 
the point where all sublattice magnetizations vanish. 

Near T c ($ M -> oo and (S*) -> 0) Eq. (jA - !^ ca n be 
simplified. Expanding in <j> M and using Eq. |AT3|) one 
obtains 

(A14) 



From Eq. I|A14(1 . it follows that for spin-independent 
Heisenberg exchange parameters [Eq. (| Alfl ] the depen- 
dence of the Curie temperature on the spin value is de- 
fined by the factor (5 M + 1)/-S>- 



The classical limit can be obtained by letting — » 
oo iriEqs. jSE} and l(AT4|4^ Factor (5, t + in 
Eq. (| A 14(1 becomes in this case unity. The temperature 
dependence of the magnetization can be calculated using 
a semiclassical analog of Eq. (| Al 2|> given by^£2 



^ ^ q.v * ' 

. . (A15) 

where C(x) — coth(x) — 1/x is the Langevin function and 
e M is the angular momentum vector of size one. 
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